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Abstract. Recent studies suggest that when cold nu- 
clear matter is compressed to high nuclear densities, di- 
quarks with spin zero and antisymmetric color wave func- 
tion Bose condensate into a superfluid/superconducting 
state that is several times as dense. Various astrophysi- 
cal phenomena may be explained by gravitational collapse 
of neutron stars (NSs) to (di)quark stars (QSs) as a re- 
sult of a first order phase transition in NSs within ~ 10^ 
years after their birth in supernova explosions, when they 
cooled and spun down sufficiently (by magnetic braking 
?). The gravitational energy release drives an explosion 
which may eject both highly relativistic narrowly coUi- 
mated jets and a mildly relativistic "spherical" shell. The 
slow contraction/cooling of the remnant QSs can power 
soft gamma ray repeaters (SGRs) and anomalous X-ray 
pulsars (AXPs), without invoking a huge magnetic energy 
storage. The jets can produce the observed gamma ray 
bursts (GRBs) in distant galaxies when they happen to 
point in our direction and the GRBs afterglows. Also the 
expanding shell from the explosion may produce a GRB 
afterglow. The jets distort the original SNR, sweep up am- 
bient matter along their trajectories, accelerate it to cos- 
mic ray (CR) energies and disperse it in hot spots which 
they form when they stop in the galactic halo. Such events 
in our Galaxy may be the main source of Galactic cosmic 
rays at all energies. 



1. First Order Phase Transition in CooHng NSs ? 

The properties of hadronic matter under extreme condi- 
tions is subject to intense theoretical studies and experi- 
mental investigations in high energy heavy ion collisions. 
Numerical techniques cannot yet simulate the behavior of 
quarks at very high densities, but, theoretical studies (Al- 
ford et al. 1998; Rapp et al. 1998; Wilczek 1998) suggest 
that in cold (T < Tc ^ 50 MeV) and dense nuclear matter 
(n > 1 fm~^) up and down quarks produce Cooper pairs 
with spin zero and antisymmetric color wave function 
which Bose condensate into a superconducting/superfluid 
ground state. The heat release in energetic heavy ion colli- 
sions limits the experimental information on dense nuclear 



matter to very high temperatures. Information on the be- 
havior of cold nuclear matter at very high densities must 
be extracted from observations of neutron stars. Here, I 
outline a possible connection between a first order phase 
transition of nuclear matter to diquark nratter in NSs, 
GRBs, SGRs and galactic CRs. 

2. Soft Gamma Ray Repeaters 

SGRs are slowly rotating (P ^^ 5 — 8 s) neutron stars that 
produce multiple bursts of soft 7-rays, often at super- 
Eddington luminosities (see, e.g.. Hurley and Kouveliotou 
in these proceedings). Four of these objects have been dis- 
covered. Two of these objects (SGR 0526-66 and SGR 
1900-1-14) have also produced hard and extremely intense 
bursts (lO^'*-^ erg in a few tenths of a second, assuming 
isotropic emission). The period, P=5.16 s, and the pe- 
riod derivative, P = 1.1 x lO"!", of SGR 1900-hl4 (see, 
e.g., Hurley et al. 1998c and references therein) yield 
a characteristic age, r = P/2P « 750 y, much younger 
than the age of its nearby supernova remnant (SNR) 
G42.8-F0.6. If this characteristic age of 1900-^14 repre- 
sents its true age and if it was born at the center of 
the SNR it must have moved with perpendicular veloc- 
ity vj^ ^ 40000(D/5.7 kpc) km s^^ to its present loca- 
tion (Hurley et al. 1998a; 1998b). The period, P=7.47 s, 
and the period derivative, P = 8 x lO"!", of SGR 1806- 
20 (Kouveliotou et al. 1998), yield a characteristic age, 
r — P/2P « 1500 y, much younger than the age of the 
nearby supernova remnant (SNR) GlO.0-0.3. If this age 
of 1806-1-20 represents its true age and if it was born at 
the center of the SNR, it must have moved with a per- 
pendicular velocity v^ « 5500(D/14.5 kpc) km s~^, to its 
present location (Kulkarni 1994). If the characteristic age 
of SGR 0526-66 in SNR N49 in the LMC, which has a 
period of P ^ 8 s, is similar to that of SGR 1900-1-14, i.e., 
T ^ 750 y, than it is also much younger than the age of 
SNR N49, and if it was born at the center of the N49 in 
the LMC, it must have moved with a perpendicular ve- 
locity vj^ r^ 5000 — 20000 km s^^ to its present location 
(Cline et al. 1982). These velocities are larger than those 
of young pulsars, such as the Crab and Vela pulsars, and 
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those of millisecond pulsars (e.g., Toscano et al 1998) by 
about two orders of magnitude. A similar situation exists 
for anomalous X-ray pulsars (AXPs). 

The inferred ages and velocities of SGRs and AXPs 
suggest that their characteristic spin down ages are not 
their true ages. I propose that they are the characteris- 
tic times that these pulsars spent in their present phase, 
a young (di)quark star (QS) that was formed recently, 
together with a GRB, by a first order phase transi- 
tion in neutron stars (NSs) that have cooled and spun 
down sufficiently. Their typical natal kick velocities are 
V ~ 500 km s~^, much larger than the typical velocities, 
(^ 100 km s~^), of young neutron stars produced in SNe, 
such as the Crab (Caraveo and Mignani 1999) and the 
Vela (Nasuti et al. 1997) pulsars and of all millisecond 
pulsars (Toscano et al. 1998). Because of their high natal 
velocities, they are expected to be found alone or in wide 
binaries, but not in close binaries. Their slow rotation may 
result from concentration of their angular momentum in 
internal superfluid vortices which are pinned to their crust 
and brake its rotation. The quiescent thermal X-ray emis- 
sion from the cooling QS may be used to verify that their 
radii are significantly smaller than those of millisecond 
pulsars. Finally, the QS end up as slowly rotating normal 
radio pulsars. 

3. Highly Relativistic Jets Prom NS Collapse 

Relativistic jets seem to be emitted when mass is ac- 
creted at a high rate onto compact objects. They are ob- 
served in galactic superluminal sources, such as the mi- 
croquasar GRS 1915-1-105 (Rodriguez and Mirabel 1994), 
where mass is accreted onto a stellar black hole (BH), in 
star binaries containing NSs such as SS433 (e.g., Hjellm- 
ing & Johnston 1988), and in many active galactic nuclei 
(AGN), where mass is accreted onto a supermassive BH. 
The ejection of the relativistic jets is not well understood 
yet. But, it seems very likely that highly relativistic jets 
are also ejected in NSs collapse because the accretion rates 
and the magnetic fields involved in such collapse are much 
larger. We estimate their kinetic energy Ejot and bulk mo- 
tion Lorentz factor F as follows: 

Most of the gravitational energy release in the gravi- 
tational collapse of NS (M ~ I.4M0) to QS, is radiated 
in a short burst of neutrinos and only a small part of it 
is released as kinetic energy of the remnant QS and of 
the ejecta which probably consists of a mildly relativistic 
spherical shell and highly relativistic two antiparralel jets. 
If momentum imbalance in the ejection of the relativistic 
jets ( and not asymmetric neutrino emission) is responsible 
for the observed large mean velocity (Lyne and Lorimcr 
1994) V « 450 ± 90 km s~^, of slowly spinning pulsars 
(presumably DQ), then momentum conservation implies 
that the kinetic energy of the jets satisfies 



If two antiparallel jets are ejected, then the jet kinetic 
energy may be of the order of Ejot ~ 10^^ erg. 

High-resolution radio observations have resolved nar- 
rowly coUimated relativistic jets from microquasars, 
quasars and radio galaxies into clouds of plasma (plas- 
moids) that are emitted in injection episodes which are 
correlated with sudden removal of the accretion disk ma- 
terial. In GRS 1915+105 these plasmoids appear to ex- 
pand freely in the plasmoid rest frame, (i.e., with the 
speed of sound in a relativistic gas, Cg ^ c/\/3, corre- 
sponding to vertical expansion speed of ^ c/Fv3 in the 
lab frame), up to a radius (Rodriguez and Mirabel 1998) 
Rp ~ 2 X 10^^ cm and retain a constant or slowly expand- 
ing radius afterwards, probably, due to magnetic con- 
finement. Likewise plasmoids from quasars seem to re- 
tain a constant radius after initial expansion. The highly 
relativistic magnetized plasmoids from NSs collapse are 
strong emitters of beamed 7-rays through synchrotron 
emission, inverse Compton scattering and resonance scat- 
tering of interstellar light. When they point in our di- 
rection from external galaxies, they can produce the ob- 
served GRBs and their afterglows (see next section). If 
the true rate of GRBs is similar to the birth rate of NSs, 
Nns ~ 2 X 10^^ y^^ in galaxies like our own (e.g., van den 
Bergh and Tamman 1991), then, since the estimated cur- 
rent rate of GRBs in galaxies like our own is (e.g., Wijers 
et al. 1997) Ngrb ^ 10^* y^^, it implies that the radiation 
emitted in GRBs must be narrowly beamed into a solid 
angle Ail that satisfies 



An ~ 27rNgrb/Nn 



TT X 10" 



(2) 



Ejet > cP„, - vM„,c - 4 X 10-^^ erg 



(1) 



GRB observations and the recent redshift measure- 
ments/estimates of some GRBs and of host galaxies of 
GRBs (see this volume) do suggest that GRBs are pro- 
duced by highly relativistic, narrowly collimated jets (Sha- 
viv and Dar 1995; Dar 1998 and references therein) with 
bulk motion velocities v corresponding to Lorentz factors 
F ^ 10"^ (Baring and Harding 1997). Emission from nar- 
row jets with F = 10'^, indeed, is beamed into Afi ~ 
tt/F^ ~ tt X 10^® consistent with eq. 2. Such strong beam- 
ing implies that we observe only a very small fraction of 
a large rate of the events that can produce GRBs. We 
can these events "Galactic" GRBs (GGRBs) if they oc- 
cur in our Milky Way (MW) galaxy and Cosmological 
GRBs (GGRBs) if they occur in distant galaxies. Eqs. 
1,2 also imply that the ejected jet (plasmoid) has a mass 

Mjet ~ 1.5 X 10-6M„, ~ 2.1 X IO-6M0 ~ 0.7MEarth. 

Even if only a fraction 77 ^ 10^^ of the jet kinetic energy is 
radiated in 7-rays, the inferred "isotropic" 7-ray emission 
in GRBs is Eisot — 'iv^^^ict ^ 4 x 10"''^ erg, while the true 
7-ray emission is only E^ ^ 10*^ erg. 



4. Jetted GRBs 

Observations of GRB afterglows suggest that GRBs are 
produced in star burst regions (SBR) where the column 
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densities of gas and light are high. The ejected jets consist 
of pure e"'"e~ plasmoids or of normal hadronic gas/plasma 
clouds. The GRB can be produced by electron synchrotron 
emission (for various GRB scenarios based on synchrotron 
emission see papers in this volume and references therein) . 
If the jet consist of a single plasmoid, then individual 7- 
ray pulses that make GRBs may be produced by magnetic 
shocks induced by internal instabilities or by external in- 
homegeneities. If the jet consists of multiple ejections of 
plasmoids, then the GRB pulses may be produced when 
later ejected plasmoids collide with earlier ejected plas- 
moids that have been slowed down by sweeping up the 
interstellar medium in front of them. Perhaps, Early time 
variability of GRBs afterglows may be used to distinguish 
between the single and multiple ejection scenarios. How- 
ever, these scenarios do not seem to provide a simple ex- 
planation for the ~ MeV "peak energy" of GRBs. Other 
emission mechanisms do provide such an explanation: 

If the highly relativistic plasmoid consists of a pure 
e"'"e~ plasma, then inverse Compton scattering of stellar 
light (hz/ = eov x 1 cV) by the plasmoid can explain the 
observed typical 7 energy (e-y ^ 4r|eov/3(l -I- z) MeV), 
GRB duration (T ~ RsBR/2cr^ ^50 s), pulse duration 
(tp - Rp/2cr2 - 150 ms), fluence (F^, - 10"^ erg cm-^), 
light curve and spectral evolution of GRBs (Shaviv and 
Dar 1995; Shaviv 1996; Dar 1998). For instance. 



o-TNeT,Ejet(l + z) 10-''^Z2N22r3e;E52 erg 



rmpc2 B'^An 



D^ 



where D = D29 x 10^^ cm is the luminosity distance of the 
GRB at redshift z, Z2 = (1 + z)/2, N = N22 x 10^^ cm^^ 
is the column density of photons along the jet trajectory in 
the star burst region, a.^, = 0.65 x 10"^'* cm^ is the Thom- 
son cross section, Ejet = E52 x 10^^ erg and F = Fs x 10^. 
If the plasmoid consists of normal crust material of 
neutron stars, (Doppler shifted Kq iron line was detected 
from the jets of SS433), than photoabsorption of stellar 
light by partially ionized iron (Doppler shifted to X-rays in 
the jet rest frame) and its reemission as 7 rays (iron X-rays 
lines in the jet rest frame) yield e-y = Fex/(1 + z) ^ McV) 
in the observer frame and 



may also produce additional (stationary) afterglow (many 
planetary nebulae, e.g., NGC 7009, NGC 6826, and some 
SNR look like ejection of antiparralel jets from a spher- 
ical explosion). When transformed to the rest frame of 
the decelerating plasmoid, the synchrotron spectra can be 
modeled by convolving typical electron energy spectrum 
(E~P at low energies up to some "break energy" where 
it steepens to E^p^^ and cutoffs exponentially at some 
higher energy due to synchrotron losses in Magnetic accel- 
eration) with the synchrotron Green's function (see, e.g., 
Meisenheimer et al. 1989). In the observer frame it yields 
spectral intensity (Dar 1998) 



I. 



^."t^ 



-0.75±0.25.-1.25±0.08 



(5) 



where a = — (p — l)/2 and (3 — — (p + 5)/6 and where I 
assumed p=: 2.5 ±0.5 for Magnetic Fermi acceleration. 
This prediction is in agreement with observations of GRB 
afterglows. Moreover, the glows of microquasar plasmoids 
and radio quasar jets after ejection and of blazar jets af- 
ter flares show the same universal behaviour as observed 
in GRBs afterglows. For instance, the glows of the ejected 
plasmoids from GRS 1915+105 on April 16, 1994 near the 
source had a — —0.8 ±0.1 and f3 = —1.3 ±0.2 (Rodriguez 
and Mirabel 1998) identical to those observed for SS 433 
(HjcUming Johnston 1988) and for the inner regions of 
jets of some radio galaxies (Bridle & Perley 1984). 



(3) 6. GGRBs - The Source of Cosmic Rays 



CTaNe^ Ejet(l 



10 ^Z2Cri9N22exr3E52 crg 



^ ~ FMFeC2 



D2A1] 



D2 

^29 



cm^ 



(4) 



where a^ = cig x 10 ^^ cm^ is the mean photoabsorption 
cross section of X-rays by partially ionized iron. 

5. GRB Afterglov^fs 

The afterglows of GRBs may be synchrotron emission 
from the decelerating plasmoids (e.g., Chiang and Der- 
mer 1997), and then they are highly beamed and may ex- 
hibit detectable superluminal velocities (v^ < Fc during 
and right after the GRB). The deceleration of a mildly 
relativistic spherical ejecta from the NS collapse to QS, 



According to the current paradigm of CR origin, CR nuclei 
with energies below 3 • 10"'^^ eV (the "knee") are acceler- 
ated in Galactic SNRs (Ginzburg 1957) and those above 
3 X 10^* eV (the "ankle" ), for which a disk origin is unlikely 
due to their isotropy, in sources far beyond our Galaxy 
(Burbidge 1964). However, recent observations, in partic- 
ular the absence (Takeda 1998) of the "GZK cutoff" in the 
CR intensity at energies above ~ lO^" eV due to interac- 
tions with the cosmic microwave radiation (Greisen 1996; 
Zatsepin & Kuz'min 1966) have put into question (e.g., 
Hillas 1998) the current paradigm for the CR origin. Dar 
and Plaga (1998) have recently proposed an alternative 
paradigm where Galactic GRBs (GGRBs) are the main 
source of the CRs at all energies and no GZK cutoff is 
expected: 

Relativistic jets are efficient CR accelerators (e.g., 
Dar 1998b). The highly relativistic narrowly coUimated 
jets/plasmoids from the birth or collapse of NSs in the 
disk of our Galaxy that are emitted with Ejot ^ 10^^ erg 
perpendicular to the Galactic disk, stop only in the 
Galactic halo, when the rest mass energy of the swept 
up ambient material becomes > their initial kinetic en- 
ergy. Through the Fermi mechanism, they accelerate the 
swept up ambient matter to CR energies and disperse 
it into the halo from the hot spots which they form 
when they finally stop in the Galactic halo (Fig.l). The 
typical equipartition magnetic fields in such hot spots 
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may reach B ~ (3Ejct/Rp)^^^ ~ 1 G. Synchrotron losses 
cut off Fermi acceleration of CR nuclei with mass num- 
ber A at E-rA2z-3/2(B/G)"i/2 X 1020 eV. Particle- 
escape cuts off Fermi acceleration when the Larmor ra- 
dius of the accelerated particles in the plasmoid rest 
frame becomes comparable to the radius of the plasmoid, 
i.e., above E ~ rZ(B/G)(Rp/0.1 pc) x lO^^ eV. Conse- 
quently, CR with E > Z X 10^" eV can no longer be 
isotropized by acceleration or deflection in hot spots with 
r ~ 1. Thus, according to Dar and Plaga (1998) CR with 
energies above lO^" eV, probably, are heavy nuclei. 

Fermi acceleration in or by the highly relativistic 
jets from GRBs can produce a broken power-law spec- 
trum, dn/dE ~ E^'^, with (3 ~ 2.2 below a knee around 
Ekncc ^ A PcV and /3 ^ 2.5 above this energy (Dar 
1998b). Spectral indices f3 ~ 2.2 were obtained also in nu- 
merical simulation of relativistic shock acceleration (e.g., 
Bednarz and Ostrowski 1998) . Galactic magnetic con- 
finement increases the density of Galactic CR by the ratio 
CTh/Rc where Th(E), is the mean residence time in the 
halo of Galactic CR with energy E and Rq ~ 50 kpc is 
the radius of the Galactic magnetic-confinement region. 
With the standard choice for the energy dependence of 
the diffusion constant (observed e.g., in solar-system plas- 
mas) one gets: Th ex (E/Z)~°'^. Consequently, the energy 
spectrum of CR is predicted to be 



dn/dE ^ C(E/Ekncc)^" 



(6) 



with a ~ (3 + 0.5 ~ 2.7 (~ 3) below (above) the knee. This 
power-law continues as long as the Galactic magnetic field 
confines the CRs. 

Part of the kinetic energy released by GGRBs is 
transported into the Galactic halo by the jets. Assuming 
equipartition of this energy, without large losses, between 
CR, gas and magnetic fields in the halo during the resi- 
dence time of CR there, the magnetic field strength Bjj in 
the halo is expected to be comparable to that of the disk 
Bh ~ (2LMw[CR]Th/R3)i/2 ~ 3 AiG where Th - 5 x 10^ 
is the mean residence time of the bulk of the CR in the 
Galactic halo. Cosmic rays with Larmor radius larger than 
the coherence length A of the halo magnetic fields, i.e., 
with energy above 



E, 



3 X 10^«(ZBh/3/iG)(A/kpc) eV, 



(7) 



escape Galactic trapping. Thus, the CR ankle is explained 
as the energy where the mean residence time Th(E) of CR 
becomes comparable to the free escape time from the halo 
Tfreo "^ 1.6(i?h/50 kpc) X 10^ ycars. Therefore, the spec- 
trum of CR with energies above the ankle, that do not 
suffer Galactic magnetic trapping, is the CR spectrum 
produced by the jet, i.e., 

dn/dE ~ C(Eanklc/Ekncc)" (E/Eanklc) ' ; E > Eanklc-(8) 

Eqs. 6,7,8 describe well the overwhole CR energy spec- 
trum. 




Fig. 1. A highly schematic sketch of the Dar-Plaga 
paradigm. The birth of QS or the collapse of NSs in the 
disk of our Galaxy leads to an ejection of two opposite jets 
that produce "hot spots" when they stop in an extended 
Galactic halo. 
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